Description of the Data The dataset contains detailed information on reported crimes, categorized by type, location, and time. It includes columns such as:
Date: The date of the occurrence. Time: The time of the occurrence. Latitude and Longitude: Geographical coordinates of the incident. Offense Category: The type of crime, such as larceny, vandalism, or burglary. Violent/Property Crime Indicator: Classification of crimes as violent or property-related. Domestic: Indicates if the crime was domestic in nature. Total Incidents: Number of incidents reported for a specific category. The data spans multiple years and includes seasonality, trend, and random components. By analyzing these elements, we aim to uncover patterns and develop reliable forecasting models to predict future crime occurrences.
Objective: Using advanced time-series methods such as ARIMA, STL decomposition, and Holt-Winters, the analysis aims to identify the most accurate model for predicting crime trends in various categories, thus enabling actionable insights.Importance: Accurately forecasting crime trends is vital for law enforcement, policymakers, and communities. By understanding the future trajectory of specific crime categories, stakeholders can proactively address challenges, optimize resource allocation, and implement preventive measures. This analysis supports data-driven decision-making, enhancing public safety and fostering community well-being.
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
##
## Attaching package: 'lubridate'
## The following objects are masked from 'package:base':
##
## date, intersect, setdiff, union
## Registered S3 method overwritten by 'quantmod':
## method from
## as.zoo.data.frame zoo
## Warning: package 'Metrics' was built under R version 4.4.2
##
## Attaching package: 'Metrics'
## The following object is masked from 'package:forecast':
##
## accuracy
We would cut the data at 2001
New
column names to be applied Date Time Lat Lon Cat VP Domestic Tot_Inc
19118 1999-12-31 19:30:00 NA NA Larceny Part I N 1 37292 1999-12-31
20:30:00 NA NA Larceny Part I N 1 119123 1999-12-31 20:00:00 NA NA
Vandalism Part II N 1 9379 2000-01-01 12:01:00 NA NA Larceny Part I N 1
10733 2000-01-01 11:00:00 NA NA Larceny Part I N 1 20324 2000-01-01
12:01:00 NA NA Forgery Part II N 1
density
plot for crime rate based on lattitude and longitude
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Warning: The dot-dot notation (`..level..`) was deprecated in ggplot2 3.4.0.
## ℹ Please use `after_stat(level)` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Warning: Removed 1 row containing non-finite outside the scale range
## (`stat_align()`).
## Warning: Removed 1 row containing missing values or values outside the scale range
## (`geom_line()`).
There appears to be seasonality in the time series, indicating the
increase of crime in specific
days.
The
residuals are randomly distributed, shows no pattern. This means that
the Naive forecast is being able to capture the seasonality and trend
Since
the histogram is normal, it means that the Naive Forecast might be a
good method in this case.
Forecast
assumes future values equal the mean of past values.
Observing
the residuals, we can say the Mean Forecast doesnot effectively captures
the pattern. Since the existence of patterns in residuals, it suggest
the models limitations.
Since the residuals do not show patterns, this suggests that the model
is capturing the seasonality or trend.
Similar
to the one without drift , residuals are randomly distributed; absence
of patterns suggest no accounted trends or seasonality.
SES
can be effective because it uses a weighted average of past values,
putting more weight on recent observations without explicitly modeling
trend or seasonality.
Since there is a significant overlap between the original timeseries and
seasonally independent timeseries, seasonality is not very pronounced.
Factors like trend or randomness moght have bigger influence on the
data.
[1]
“Number of differences required for stationarity: 1”
Fitting models using approximations to speed things up…
ARIMA(0,1,0) : 3247.037 ARIMA(0,1,0) with drift : 3249.007 ARIMA(0,1,0)(0,0,1)[12] : 3237.545 ARIMA(0,1,0)(0,0,1)[12] with drift : 3239.523 ARIMA(0,1,0)(0,0,2)[12] : 3237.29 ARIMA(0,1,0)(0,0,2)[12] with drift : 3239.277 ARIMA(0,1,0)(1,0,0)[12] : 3204.453 ARIMA(0,1,0)(1,0,0)[12] with drift : 3206.394 ARIMA(0,1,0)(1,0,1)[12] : 3204.482 ARIMA(0,1,0)(1,0,1)[12] with drift : 3206.477 ARIMA(0,1,0)(1,0,2)[12] : 3203.137 ARIMA(0,1,0)(1,0,2)[12] with drift : 3205.167 ARIMA(0,1,0)(2,0,0)[12] : 3204.508 ARIMA(0,1,0)(2,0,0)[12] with drift : 3206.511 ARIMA(0,1,0)(2,0,1)[12] : Inf ARIMA(0,1,0)(2,0,1)[12] with drift : Inf ARIMA(0,1,0)(2,0,2)[12] : Inf ARIMA(0,1,0)(2,0,2)[12] with drift : Inf ARIMA(0,1,1) : 3199.389 ARIMA(0,1,1) with drift : 3201.038 ARIMA(0,1,1)(0,0,1)[12] : 3190.258 ARIMA(0,1,1)(0,0,1)[12] with drift : 3191.986 ARIMA(0,1,1)(0,0,2)[12] : 3188.588 ARIMA(0,1,1)(0,0,2)[12] with drift : 3190.33 ARIMA(0,1,1)(1,0,0)[12] : 3125.141 ARIMA(0,1,1)(1,0,0)[12] with drift : 3127.194 ARIMA(0,1,1)(1,0,1)[12] : 3121.32 ARIMA(0,1,1)(1,0,1)[12] with drift : 3123.372 ARIMA(0,1,1)(1,0,2)[12] : 3111.57 ARIMA(0,1,1)(1,0,2)[12] with drift : 3113.577 ARIMA(0,1,1)(2,0,0)[12] : 3117.113 ARIMA(0,1,1)(2,0,0)[12] with drift : 3119.146 ARIMA(0,1,1)(2,0,1)[12] : Inf ARIMA(0,1,1)(2,0,1)[12] with drift : Inf ARIMA(0,1,1)(2,0,2)[12] : Inf ARIMA(0,1,1)(2,0,2)[12] with drift : Inf ARIMA(0,1,2) : 3201.152 ARIMA(0,1,2) with drift : 3202.786 ARIMA(0,1,2)(0,0,1)[12] : 3192.225 ARIMA(0,1,2)(0,0,1)[12] with drift : 3193.958 ARIMA(0,1,2)(0,0,2)[12] : 3190.614 ARIMA(0,1,2)(0,0,2)[12] with drift : 3192.364 ARIMA(0,1,2)(1,0,0)[12] : 3125.941 ARIMA(0,1,2)(1,0,0)[12] with drift : 3128.009 ARIMA(0,1,2)(1,0,1)[12] : 3121.431 ARIMA(0,1,2)(1,0,1)[12] with drift : 3123.428 ARIMA(0,1,2)(1,0,2)[12] : 3112.274 ARIMA(0,1,2)(1,0,2)[12] with drift : 3114.222 ARIMA(0,1,2)(2,0,0)[12] : 3118.248 ARIMA(0,1,2)(2,0,0)[12] with drift : 3120.265 ARIMA(0,1,2)(2,0,1)[12] : Inf ARIMA(0,1,2)(2,0,1)[12] with drift : Inf ARIMA(0,1,3) : 3202.036 ARIMA(0,1,3) with drift : 3203.624 ARIMA(0,1,3)(0,0,1)[12] : 3194.019 ARIMA(0,1,3)(0,0,1)[12] with drift : 3195.746 ARIMA(0,1,3)(0,0,2)[12] : 3192.594 ARIMA(0,1,3)(0,0,2)[12] with drift : 3194.349 ARIMA(0,1,3)(1,0,0)[12] : 3123.409 ARIMA(0,1,3)(1,0,0)[12] with drift : 3125.471 ARIMA(0,1,3)(1,0,1)[12] : 3116.114 ARIMA(0,1,3)(1,0,1)[12] with drift : 3117.889 ARIMA(0,1,3)(2,0,0)[12] : 3115.711 ARIMA(0,1,3)(2,0,0)[12] with drift : 3117.591 ARIMA(0,1,4) : 3203.986 ARIMA(0,1,4) with drift : 3205.568 ARIMA(0,1,4)(0,0,1)[12] : 3196.102 ARIMA(0,1,4)(0,0,1)[12] with drift : 3197.843 ARIMA(0,1,4)(1,0,0)[12] : 3123.281 ARIMA(0,1,4)(1,0,0)[12] with drift : 3125.352 ARIMA(0,1,5) : 3205.652 ARIMA(0,1,5) with drift : 3207.278 ARIMA(1,1,0) : 3170.336 ARIMA(1,1,0) with drift : 3172.345 ARIMA(1,1,0)(0,0,1)[12] : 3154.501 ARIMA(1,1,0)(0,0,1)[12] with drift : 3156.529 ARIMA(1,1,0)(0,0,2)[12] : 3153.158 ARIMA(1,1,0)(0,0,2)[12] with drift : 3155.195 ARIMA(1,1,0)(1,0,0)[12] : 3142.964 ARIMA(1,1,0)(1,0,0)[12] with drift : 3144.987 ARIMA(1,1,0)(1,0,1)[12] : 3144.872 ARIMA(1,1,0)(1,0,1)[12] with drift : 3146.915 ARIMA(1,1,0)(1,0,2)[12] : 3142.916 ARIMA(1,1,0)(1,0,2)[12] with drift : 3144.992 ARIMA(1,1,0)(2,0,0)[12] : 3145.232 ARIMA(1,1,0)(2,0,0)[12] with drift : 3147.294 ARIMA(1,1,0)(2,0,1)[12] : Inf ARIMA(1,1,0)(2,0,1)[12] with drift : Inf ARIMA(1,1,0)(2,0,2)[12] : Inf ARIMA(1,1,0)(2,0,2)[12] with drift : Inf ARIMA(1,1,1) : 3128.264 ARIMA(1,1,1) with drift : 3130.011 ARIMA(1,1,1)(0,0,1)[12] : 3116.125 ARIMA(1,1,1)(0,0,1)[12] with drift : 3117.998 ARIMA(1,1,1)(0,0,2)[12] : 3109.291 ARIMA(1,1,1)(0,0,2)[12] with drift : 3111.139 ARIMA(1,1,1)(1,0,0)[12] : 3111.874 ARIMA(1,1,1)(1,0,0)[12] with drift : 3113.61 ARIMA(1,1,1)(1,0,1)[12] : Inf ARIMA(1,1,1)(1,0,1)[12] with drift : Inf ARIMA(1,1,1)(1,0,2)[12] : Inf ARIMA(1,1,1)(1,0,2)[12] with drift : Inf ARIMA(1,1,1)(2,0,0)[12] : 3113.181 ARIMA(1,1,1)(2,0,0)[12] with drift : 3114.944 ARIMA(1,1,1)(2,0,1)[12] : Inf ARIMA(1,1,1)(2,0,1)[12] with drift : Inf ARIMA(1,1,2) : 3129.316 ARIMA(1,1,2) with drift : 3130.953 ARIMA(1,1,2)(0,0,1)[12] : 3117.513 ARIMA(1,1,2)(0,0,1)[12] with drift : 3119.313 ARIMA(1,1,2)(0,0,2)[12] : 3111.201 ARIMA(1,1,2)(0,0,2)[12] with drift : 3113.025 ARIMA(1,1,2)(1,0,0)[12] : 3113.703 ARIMA(1,1,2)(1,0,0)[12] with drift : 3115.39 ARIMA(1,1,2)(1,0,1)[12] : Inf ARIMA(1,1,2)(1,0,1)[12] with drift : Inf ARIMA(1,1,2)(2,0,0)[12] : 3114.993 ARIMA(1,1,2)(2,0,0)[12] with drift : 3116.706 ARIMA(1,1,3) : 3125.364 ARIMA(1,1,3) with drift : 3126.852 ARIMA(1,1,3)(0,0,1)[12] : 3112.987 ARIMA(1,1,3)(0,0,1)[12] with drift : 3114.59 ARIMA(1,1,3)(1,0,0)[12] : 3109.521 ARIMA(1,1,3)(1,0,0)[12] with drift : 3110.98 ARIMA(1,1,4) : 3126.859 ARIMA(1,1,4) with drift : 3128.277 ARIMA(2,1,0) : 3156.71 ARIMA(2,1,0) with drift : 3158.725 ARIMA(2,1,0)(0,0,1)[12] : 3136.907 ARIMA(2,1,0)(0,0,1)[12] with drift : 3138.943 ARIMA(2,1,0)(0,0,2)[12] : 3132.097 ARIMA(2,1,0)(0,0,2)[12] with drift : 3134.141 ARIMA(2,1,0)(1,0,0)[12] : 3126.481 ARIMA(2,1,0)(1,0,0)[12] with drift : 3128.549 ARIMA(2,1,0)(1,0,1)[12] : Inf ARIMA(2,1,0)(1,0,1)[12] with drift : Inf ARIMA(2,1,0)(1,0,2)[12] : Inf ARIMA(2,1,0)(1,0,2)[12] with drift : Inf ARIMA(2,1,0)(2,0,0)[12] : 3129.669 ARIMA(2,1,0)(2,0,0)[12] with drift : 3131.74 ARIMA(2,1,0)(2,0,1)[12] : Inf ARIMA(2,1,0)(2,0,1)[12] with drift : Inf ARIMA(2,1,1) : 3123.856 ARIMA(2,1,1) with drift : 3125.479 ARIMA(2,1,1)(0,0,1)[12] : 3109.51 ARIMA(2,1,1)(0,0,1)[12] with drift : 3111.24 ARIMA(2,1,1)(0,0,2)[12] : 3105.398 ARIMA(2,1,1)(0,0,2)[12] with drift : 3107.135 ARIMA(2,1,1)(1,0,0)[12] : 3111.809 ARIMA(2,1,1)(1,0,0)[12] with drift : 3113.19 ARIMA(2,1,1)(1,0,1)[12] : Inf ARIMA(2,1,1)(1,0,1)[12] with drift : Inf ARIMA(2,1,1)(2,0,0)[12] : 3118.412 ARIMA(2,1,1)(2,0,0)[12] with drift : 3120.142 ARIMA(2,1,2) : 3125.515 ARIMA(2,1,2) with drift : 3127.079 ARIMA(2,1,2)(0,0,1)[12] : Inf ARIMA(2,1,2)(0,0,1)[12] with drift : Inf ARIMA(2,1,2)(1,0,0)[12] : 3113.878 ARIMA(2,1,2)(1,0,0)[12] with drift : 3115.281 ARIMA(2,1,3) : 3127.311 ARIMA(2,1,3) with drift : 3128.82 ARIMA(3,1,0) : 3154.761 ARIMA(3,1,0) with drift : 3156.761 ARIMA(3,1,0)(0,0,1)[12] : 3135.46 ARIMA(3,1,0)(0,0,1)[12] with drift : 3137.482 ARIMA(3,1,0)(0,0,2)[12] : 3129.146 ARIMA(3,1,0)(0,0,2)[12] with drift : 3131.169 ARIMA(3,1,0)(1,0,0)[12] : 3125.889 ARIMA(3,1,0)(1,0,0)[12] with drift : 3127.97 ARIMA(3,1,0)(1,0,1)[12] : Inf ARIMA(3,1,0)(1,0,1)[12] with drift : Inf ARIMA(3,1,0)(2,0,0)[12] : 3127.238 ARIMA(3,1,0)(2,0,0)[12] with drift : 3129.302 ARIMA(3,1,1) : 3128.525 ARIMA(3,1,1) with drift : 3129.79 ARIMA(3,1,1)(0,0,1)[12] : 3109.761 ARIMA(3,1,1)(0,0,1)[12] with drift : Inf ARIMA(3,1,1)(1,0,0)[12] : Inf ARIMA(3,1,1)(1,0,0)[12] with drift : Inf ARIMA(3,1,2) : 3128.865 ARIMA(3,1,2) with drift : 3130.242 ARIMA(4,1,0) : 3150.852 ARIMA(4,1,0) with drift : 3152.825 ARIMA(4,1,0)(0,0,1)[12] : 3134.169 ARIMA(4,1,0)(0,0,1)[12] with drift : 3136.199 ARIMA(4,1,0)(1,0,0)[12] : 3124.297 ARIMA(4,1,0)(1,0,0)[12] with drift : 3126.393 ARIMA(4,1,1) : 3129.882 ARIMA(4,1,1) with drift : 3131.312 ARIMA(5,1,0) : 3146.938 ARIMA(5,1,0) with drift : 3148.96
Now re-fitting the best model(s) without approximations…
Best model: ARIMA(2,1,1)(0,0,2)[12]
Series: ts_monthly ARIMA(2,1,1)(0,0,2)[12]
Coefficients: ar1 ar2 ma1 sma1 sma2 0.3127 0.2151 -0.9549 0.2564 0.1379 s.e. 0.0777 0.0771 0.0316 0.0745 0.0666
sigma^2 = 2291: log likelihood = -1579.44 AIC=3170.87 AICc=3171.16 BIC=3193.08
Training set error measures: ME RMSE MAE MPE MAPE MASE ACF1 Training
set 5.442304 47.38631 34.60418 0.3967119 8.848714 0.7650856 0.05655193
The ARIMA(2,1,1)(0,0,2)[12] model indicates a seasonal ARIMA with two
autoregressive terms (AR), one differencing step (I), and one moving
average term (MA), alongside two seasonal moving average terms (sMA) for
monthly data. The coefficients provide weights for these components,
with associated standard errors indicating the uncertainty of each
estimate. Key model evaluation metrics include AIC (3170.87), AICc
(3171.16), and BIC (3193.08), where lower values indicate better fit
while penalizing complexity. Training set error measures such as RMSE
(47.39) and MAE (34.60) highlight the average prediction errors. MAPE
(8.85%) reflects a reasonably accurate model, with residuals showing
minimal autocorrelation (ACF1 = 0.056), suggesting a good fit for the
data with little remaining structure in the errors.
The
randomness in the residuals suggest a good fit.[1] “Ljung-Box Test
Results:”
Box-Ljung test
data: ts_monthly_arima_residuals X-squared = 20.077, df = 20, p-value = 0.4531
The residuals do not show significant autocorrelation, indicating
that the ARIMA model has captured the dependencies in the data well, and
no major patterns remain unexplained. [1] “Residual Statistics:” Metric
Value 1 Mean 5.442304 2 Variance 2223.254526 3 Standard Deviation
47.151400
Method RMSE 1 Naive 55.75517 2 SES 49.67305 3 HW 42.96673 4 MA3 28.69516
5 MA6 31.16097 6 MA9 35.25139
Processing category: Larceny
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 2641.64 ARIMA(0,1,0) with drift : 2776.66 ARIMA(1,1,0)(1,0,0)[12] with drift : 2683.655 ARIMA(0,1,1)(0,0,1)[12] with drift : 2709.827 ARIMA(0,1,0) : 2774.658 ARIMA(2,1,2)(0,0,1)[12] with drift : Inf ARIMA(2,1,2)(1,0,0)[12] with drift : 2639.555 ARIMA(2,1,2) with drift : 2649.891 ARIMA(2,1,2)(2,0,0)[12] with drift : 2639.948 ARIMA(2,1,2)(2,0,1)[12] with drift : 2637.192 ARIMA(2,1,2)(2,0,2)[12] with drift : Inf ARIMA(2,1,2)(1,0,2)[12] with drift : 2640.357 ARIMA(1,1,2)(2,0,1)[12] with drift : 2636.174 ARIMA(1,1,2)(1,0,1)[12] with drift : 2640.009 ARIMA(1,1,2)(2,0,0)[12] with drift : 2638.745 ARIMA(1,1,2)(2,0,2)[12] with drift : Inf ARIMA(1,1,2)(1,0,0)[12] with drift : 2638.197 ARIMA(1,1,2)(1,0,2)[12] with drift : 2639.592 ARIMA(0,1,2)(2,0,1)[12] with drift : 2643.392 ARIMA(1,1,1)(2,0,1)[12] with drift : 2634.087 ARIMA(1,1,1)(1,0,1)[12] with drift : 2638.027 ARIMA(1,1,1)(2,0,0)[12] with drift : 2636.886 ARIMA(1,1,1)(2,0,2)[12] with drift : Inf ARIMA(1,1,1)(1,0,0)[12] with drift : 2636.14 ARIMA(1,1,1)(1,0,2)[12] with drift : 2637.483 ARIMA(0,1,1)(2,0,1)[12] with drift : 2642.964 ARIMA(1,1,0)(2,0,1)[12] with drift : 2679.545 ARIMA(2,1,1)(2,0,1)[12] with drift : 2636.359 ARIMA(0,1,0)(2,0,1)[12] with drift : 2745.633 ARIMA(2,1,0)(2,0,1)[12] with drift : 2669.133 ARIMA(1,1,1)(2,0,1)[12] : 2632.422 ARIMA(1,1,1)(1,0,1)[12] : 2636.08 ARIMA(1,1,1)(2,0,0)[12] : 2634.995 ARIMA(1,1,1)(2,0,2)[12] : Inf ARIMA(1,1,1)(1,0,0)[12] : 2634.193 ARIMA(1,1,1)(1,0,2)[12] : 2635.513 ARIMA(0,1,1)(2,0,1)[12] : 2640.881 ARIMA(1,1,0)(2,0,1)[12] : 2677.474 ARIMA(2,1,1)(2,0,1)[12] : 2634.796 ARIMA(1,1,2)(2,0,1)[12] : 2634.476 ARIMA(0,1,0)(2,0,1)[12] : 2743.635 ARIMA(0,1,2)(2,0,1)[12] : 2641.295 ARIMA(2,1,0)(2,0,1)[12] : 2667.037 ARIMA(2,1,2)(2,0,1)[12] : 2635.575
Now re-fitting the best model(s) without approximations…
ARIMA(1,1,1)(2,0,1)[12] : Inf ARIMA(1,1,1)(2,0,1)[12] with drift : Inf ARIMA(1,1,1)(1,0,0)[12] : 2683.859
Best model: ARIMA(1,1,1)(1,0,0)[12]
Ljung-Box
Test p-value for Larceny - Naive Forecast : 1.527187e-06
Ljung-Box
Test p-value for Larceny - Mean Forecast : 0
Ljung-Box
Test p-value for Larceny - Holt-Winters Forecast : 0.1010496
Ljung-Box
Test p-value for Larceny - ARIMA Forecast : 0.1524402
Ljung-Box
Test p-value for STL Residuals: Larceny : 0.001478226 Processing
category: Vandalism
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : Inf ARIMA(0,1,0) with drift : 2332.59 ARIMA(1,1,0)(1,0,0)[12] with drift : 2261.035 ARIMA(0,1,1)(0,0,1)[12] with drift : 2232.947 ARIMA(0,1,0) : 2330.566 ARIMA(0,1,1) with drift : 2237.802 ARIMA(0,1,1)(1,0,1)[12] with drift : 2206.932 ARIMA(0,1,1)(1,0,0)[12] with drift : 2213.167 ARIMA(0,1,1)(2,0,1)[12] with drift : 2198.024 ARIMA(0,1,1)(2,0,0)[12] with drift : 2196.375 ARIMA(0,1,0)(2,0,0)[12] with drift : 2325.821 ARIMA(1,1,1)(2,0,0)[12] with drift : 2195.165 ARIMA(1,1,1)(1,0,0)[12] with drift : 2193.909 ARIMA(1,1,1) with drift : 2205.946 ARIMA(1,1,1)(1,0,1)[12] with drift : 2178.175 ARIMA(1,1,1)(0,0,1)[12] with drift : 2198.938 ARIMA(1,1,1)(2,0,1)[12] with drift : 2196.419 ARIMA(1,1,1)(1,0,2)[12] with drift : 2178.444 ARIMA(1,1,1)(0,0,2)[12] with drift : 2199.567 ARIMA(1,1,1)(2,0,2)[12] with drift : 2192.485 ARIMA(1,1,0)(1,0,1)[12] with drift : 2257.381 ARIMA(2,1,1)(1,0,1)[12] with drift : 2180.963 ARIMA(1,1,2)(1,0,1)[12] with drift : 2180.14 ARIMA(0,1,0)(1,0,1)[12] with drift : 2328.577 ARIMA(0,1,2)(1,0,1)[12] with drift : 2202.586 ARIMA(2,1,0)(1,0,1)[12] with drift : 2227.003 ARIMA(1,1,1)(1,0,1)[12] : 2178.86
Now re-fitting the best model(s) without approximations…
ARIMA(1,1,1)(1,0,1)[12] with drift : Inf ARIMA(1,1,1)(1,0,2)[12] with drift : Inf ARIMA(1,1,1)(1,0,1)[12] : Inf ARIMA(1,1,2)(1,0,1)[12] with drift : 2219.14
Best model: ARIMA(1,1,2)(1,0,1)[12] with drift
Ljung-Box
Test p-value for Vandalism - Naive Forecast : 2.482997e-08
Ljung-Box
Test p-value for Vandalism - Mean Forecast : 0
Ljung-Box
Test p-value for Vandalism - Holt-Winters Forecast : 0.9007833
Ljung-Box
Test p-value for Vandalism - ARIMA Forecast : 0.9190618
Ljung-Box
Test p-value for STL Residuals: Vandalism : 0.09399927 Processing
category: Forgery
Fitting models using approximations to speed things up…
ARIMA(2,0,2)(1,0,1)[12] with non-zero mean : 1633.015 ARIMA(0,0,0) with non-zero mean : 1801.16 ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 1625.413 ARIMA(0,0,1)(0,0,1)[12] with non-zero mean : 1687.632 ARIMA(0,0,0) with zero mean : 2039.029 ARIMA(1,0,0) with non-zero mean : 1638.37 ARIMA(1,0,0)(2,0,0)[12] with non-zero mean : 1631.362 ARIMA(1,0,0)(1,0,1)[12] with non-zero mean : 1627.392 ARIMA(1,0,0)(0,0,1)[12] with non-zero mean : 1640.343 ARIMA(1,0,0)(2,0,1)[12] with non-zero mean : 1632.356 ARIMA(0,0,0)(1,0,0)[12] with non-zero mean : 1788.799 ARIMA(2,0,0)(1,0,0)[12] with non-zero mean : 1627.932 ARIMA(1,0,1)(1,0,0)[12] with non-zero mean : 1627.368 ARIMA(0,0,1)(1,0,0)[12] with non-zero mean : 1678.182 ARIMA(2,0,1)(1,0,0)[12] with non-zero mean : 1630.024 ARIMA(1,0,0)(1,0,0)[12] with zero mean : 1655.583
Now re-fitting the best model(s) without approximations…
ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 1639.918
Best model: ARIMA(1,0,0)(1,0,0)[12] with non-zero mean
Ljung-Box
Test p-value for Forgery - Naive Forecast : 0.07279791
Ljung-Box
Test p-value for Forgery - Mean Forecast : 0
Ljung-Box
Test p-value for Forgery - Holt-Winters Forecast : 0.9515707
Ljung-Box
Test p-value for Forgery - ARIMA Forecast : 0.9967535
Ljung-Box
Test p-value for STL Residuals: Forgery : 0.01603712 Processing
category: Simple Assault
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1958.233 ARIMA(0,1,0) with drift : 2113.539 ARIMA(1,1,0)(1,0,0)[12] with drift : 2022.612 ARIMA(0,1,1)(0,0,1)[12] with drift : 1954.714 ARIMA(0,1,0) : 2111.533 ARIMA(0,1,1) with drift : 1952.66 ARIMA(0,1,1)(1,0,0)[12] with drift : 1956.031 ARIMA(0,1,1)(1,0,1)[12] with drift : 1953.576 ARIMA(1,1,1) with drift : 1955.487 ARIMA(0,1,2) with drift : 1954.518 ARIMA(1,1,0) with drift : 2020.182 ARIMA(1,1,2) with drift : 1956.568 ARIMA(0,1,1) : 1950.827 ARIMA(0,1,1)(1,0,0)[12] : 1954.109 ARIMA(0,1,1)(0,0,1)[12] : 1952.868 ARIMA(0,1,1)(1,0,1)[12] : 1951.571 ARIMA(1,1,1) : 1953.46 ARIMA(0,1,2) : 1952.679 ARIMA(1,1,0) : 2018.16 ARIMA(1,1,2) : 1954.571
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1) : 1955.386
Best model: ARIMA(0,1,1)
Ljung-Box
Test p-value for Simple Assault - Naive Forecast : 6.505907e-14
Ljung-Box
Test p-value for Simple Assault - Mean Forecast : 0
Ljung-Box
Test p-value for Simple Assault - Holt-Winters Forecast : 0.8592693
Ljung-Box
Test p-value for Simple Assault - ARIMA Forecast : 0.5248824 Ljung-Box
Test p-value for STL Residuals: Simple Assault : 0.003914289 Processing
category: Burglary
## Warning in HoltWinters(ts_category): optimization difficulties: ERROR:
## ABNORMAL_TERMINATION_IN_LNSRCH
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 2108.454 ARIMA(0,1,0) with drift : 2231.45 ARIMA(1,1,0)(1,0,0)[12] with drift : 2157.879 ARIMA(0,1,1)(0,0,1)[12] with drift : 2109.538 ARIMA(0,1,0) : 2229.433 ARIMA(2,1,2)(0,0,1)[12] with drift : 2106.862 ARIMA(2,1,2) with drift : 2109.407 ARIMA(2,1,2)(0,0,2)[12] with drift : 2100.158 ARIMA(2,1,2)(1,0,2)[12] with drift : 2088.488 ARIMA(2,1,2)(2,0,2)[12] with drift : 2101.188 ARIMA(2,1,2)(2,0,1)[12] with drift : 2100.917 ARIMA(1,1,2)(1,0,2)[12] with drift : 2104.972 ARIMA(2,1,1)(1,0,2)[12] with drift : 2089.891 ARIMA(3,1,2)(1,0,2)[12] with drift : 2104.274 ARIMA(2,1,3)(1,0,2)[12] with drift : Inf ARIMA(1,1,1)(1,0,2)[12] with drift : 2115.62 ARIMA(1,1,3)(1,0,2)[12] with drift : 2101.315 ARIMA(3,1,1)(1,0,2)[12] with drift : 2103.308 ARIMA(3,1,3)(1,0,2)[12] with drift : Inf ARIMA(2,1,2)(1,0,2)[12] : 2089.12
Now re-fitting the best model(s) without approximations…
ARIMA(2,1,2)(1,0,2)[12] with drift : 2100.016
Best model: ARIMA(2,1,2)(1,0,2)[12] with drift
Ljung-Box
Test p-value for Burglary - Naive Forecast : 8.570922e-14
Ljung-Box
Test p-value for Burglary - Mean Forecast : 0
Ljung-Box
Test p-value for Burglary - Holt-Winters Forecast : 0.01535463
Ljung-Box
Test p-value for Burglary - ARIMA Forecast : 0.762681
Ljung-Box
Test p-value for STL Residuals: Burglary : 3.810994e-06 Processing
category: Fraud
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 2113.328 ARIMA(0,1,0) with drift : 2220.713 ARIMA(1,1,0)(1,0,0)[12] with drift : 2174.969 ARIMA(0,1,1)(0,0,1)[12] with drift : 2100.989 ARIMA(0,1,0) : 2218.688 ARIMA(0,1,1) with drift : 2100.173 ARIMA(0,1,1)(1,0,0)[12] with drift : 2107.186 ARIMA(0,1,1)(1,0,1)[12] with drift : 2109.114 ARIMA(1,1,1) with drift : 2101.461 ARIMA(0,1,2) with drift : 2100.544 ARIMA(1,1,0) with drift : 2167.27 ARIMA(1,1,2) with drift : 2103.235 ARIMA(0,1,1) : 2099.778 ARIMA(0,1,1)(1,0,0)[12] : 2105.895 ARIMA(0,1,1)(0,0,1)[12] : 2100.431 ARIMA(0,1,1)(1,0,1)[12] : 2107.876 ARIMA(1,1,1) : 2101.172 ARIMA(0,1,2) : 2100.298 ARIMA(1,1,0) : 2165.273 ARIMA(1,1,2) : 2102.898
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1) : 2105.206
Best model: ARIMA(0,1,1)
Ljung-Box
Test p-value for Fraud - Naive Forecast : 3.019938e-08
Ljung-Box
Test p-value for Fraud - Mean Forecast : 0
Ljung-Box
Test p-value for Fraud - Holt-Winters Forecast : 0.786637
Ljung-Box
Test p-value for Fraud - ARIMA Forecast : 0.7715797
Ljung-Box
Test p-value for STL Residuals: Fraud : 0.02391738 Processing category:
All Other Offenses
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : Inf ARIMA(0,1,0) with drift : 2365.74 ARIMA(1,1,0)(1,0,0)[12] with drift : 2303.053 ARIMA(0,1,1)(0,0,1)[12] with drift : 2248.968 ARIMA(0,1,0) : 2363.713 ARIMA(0,1,1) with drift : 2248.692 ARIMA(0,1,1)(1,0,0)[12] with drift : 2257.24 ARIMA(0,1,1)(1,0,1)[12] with drift : 2252.381 ARIMA(1,1,1) with drift : 2252.277 ARIMA(0,1,2) with drift : 2250.517 ARIMA(1,1,0) with drift : 2294.209 ARIMA(1,1,2) with drift : 2249.89 ARIMA(0,1,1) : 2246.888 ARIMA(0,1,1)(1,0,0)[12] : 2255.337 ARIMA(0,1,1)(0,0,1)[12] : 2247.113 ARIMA(0,1,1)(1,0,1)[12] : 2250.313 ARIMA(1,1,1) : 2250.351 ARIMA(0,1,2) : 2248.72 ARIMA(1,1,0) : 2292.171 ARIMA(1,1,2) : 2248.005
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1) : 2252.451
Best model: ARIMA(0,1,1)
Ljung-Box
Test p-value for All Other Offenses - Naive Forecast : 1.627309e-11
Ljung-Box
Test p-value for All Other Offenses - Mean Forecast : 0
Ljung-Box
Test p-value for All Other Offenses - Holt-Winters Forecast : 0.6371575
Ljung-Box
Test p-value for All Other Offenses - ARIMA Forecast : 0.4050735
Ljung-Box
Test p-value for STL Residuals: All Other Offenses : 0.02077459
Processing category: Drugs
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 2036.727 ARIMA(0,1,0) with drift : 2135.05 ARIMA(1,1,0)(1,0,0)[12] with drift : 2079.031 ARIMA(0,1,1)(0,0,1)[12] with drift : 2026.478 ARIMA(0,1,0) : 2133.022 ARIMA(0,1,1) with drift : 2025.643 ARIMA(0,1,1)(1,0,0)[12] with drift : 2034.695 ARIMA(0,1,1)(1,0,1)[12] with drift : 2035.297 ARIMA(1,1,1) with drift : 2026.898 ARIMA(0,1,2) with drift : 2024.925 ARIMA(0,1,2)(1,0,0)[12] with drift : 2032.814 ARIMA(0,1,2)(0,0,1)[12] with drift : 2025.027 ARIMA(0,1,2)(1,0,1)[12] with drift : 2033.065 ARIMA(1,1,2) with drift : 2028.947 ARIMA(0,1,3) with drift : 2024.874 ARIMA(0,1,3)(1,0,0)[12] with drift : 2032.714 ARIMA(0,1,3)(0,0,1)[12] with drift : 2024.82 ARIMA(0,1,3)(1,0,1)[12] with drift : 2033.141 ARIMA(0,1,3)(0,0,2)[12] with drift : 2026.25 ARIMA(0,1,3)(1,0,2)[12] with drift : 2035.035 ARIMA(1,1,3)(0,0,1)[12] with drift : 2028.703 ARIMA(0,1,4)(0,0,1)[12] with drift : 2026.474 ARIMA(1,1,2)(0,0,1)[12] with drift : 2028.794 ARIMA(1,1,4)(0,0,1)[12] with drift : 2029.596 ARIMA(0,1,3)(0,0,1)[12] : 2022.988 ARIMA(0,1,3) : 2023.072 ARIMA(0,1,3)(1,0,1)[12] : 2031.341 ARIMA(0,1,3)(0,0,2)[12] : 2024.41 ARIMA(0,1,3)(1,0,0)[12] : 2030.959 ARIMA(0,1,3)(1,0,2)[12] : 2033.225 ARIMA(0,1,2)(0,0,1)[12] : 2023.173 ARIMA(1,1,3)(0,0,1)[12] : 2026.757 ARIMA(0,1,4)(0,0,1)[12] : 2024.642 ARIMA(1,1,2)(0,0,1)[12] : 2026.824 ARIMA(1,1,4)(0,0,1)[12] : 2027.661
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,3)(0,0,1)[12] : 2028.217
Best model: ARIMA(0,1,3)(0,0,1)[12]
Ljung-Box
Test p-value for Drugs - Naive Forecast : 5.70779e-11
Ljung-Box
Test p-value for Drugs - Mean Forecast : 0
Ljung-Box
Test p-value for Drugs - Holt-Winters Forecast : 0.05942824
Ljung-Box
Test p-value for Drugs - ARIMA Forecast : 0.3266919
Ljung-Box
Test p-value for STL Residuals: Drugs : 3.092142e-06 Processing
category: Missing Person
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1233.49 ARIMA(0,1,0) with drift : 1356.371 ARIMA(1,1,0)(1,0,0)[12] with drift : 1284.82 ARIMA(0,1,1)(0,0,1)[12] with drift : 1237.582 ARIMA(0,1,0) : 1354.361 ARIMA(2,1,2)(0,0,1)[12] with drift : Inf ARIMA(2,1,2)(1,0,0)[12] with drift : 1231.581 ARIMA(2,1,2) with drift : Inf ARIMA(2,1,2)(2,0,0)[12] with drift : Inf ARIMA(2,1,2)(2,0,1)[12] with drift : Inf ARIMA(1,1,2)(1,0,0)[12] with drift : Inf ARIMA(2,1,1)(1,0,0)[12] with drift : 1232.286 ARIMA(3,1,2)(1,0,0)[12] with drift : 1230.621 ARIMA(3,1,2) with drift : 1228.419 ARIMA(3,1,2)(0,0,1)[12] with drift : Inf ARIMA(3,1,2)(1,0,1)[12] with drift : 1232.706 ARIMA(3,1,1) with drift : 1230.742 ARIMA(4,1,2) with drift : 1224.156 ARIMA(4,1,2)(1,0,0)[12] with drift : 1232.629 ARIMA(4,1,2)(0,0,1)[12] with drift : 1225.037 ARIMA(4,1,2)(1,0,1)[12] with drift : 1231.606 ARIMA(4,1,1) with drift : 1222.09 ARIMA(4,1,1)(1,0,0)[12] with drift : 1230.498 ARIMA(4,1,1)(0,0,1)[12] with drift : 1222.934 ARIMA(4,1,1)(1,0,1)[12] with drift : 1230.422 ARIMA(4,1,0) with drift : 1230.514 ARIMA(5,1,1) with drift : 1225.051 ARIMA(3,1,0) with drift : 1253.004 ARIMA(5,1,0) with drift : 1223.072 ARIMA(5,1,2) with drift : Inf ARIMA(4,1,1) : 1219.997 ARIMA(4,1,1)(1,0,0)[12] : 1228.402 ARIMA(4,1,1)(0,0,1)[12] : 1220.824 ARIMA(4,1,1)(1,0,1)[12] : 1228.289 ARIMA(3,1,1) : 1228.767 ARIMA(4,1,0) : 1228.44 ARIMA(5,1,1) : 1222.949 ARIMA(4,1,2) : 1222.049 ARIMA(3,1,0) : 1250.958 ARIMA(3,1,2) : 1226.467 ARIMA(5,1,0) : 1220.987 ARIMA(5,1,2) : Inf
Now re-fitting the best model(s) without approximations…
ARIMA(4,1,1) : 1230.212
Best model: ARIMA(4,1,1)
Ljung-Box
Test p-value for Missing Person - Naive Forecast : 1.607603e-12
Ljung-Box
Test p-value for Missing Person - Mean Forecast : 0
Ljung-Box
Test p-value for Missing Person - Holt-Winters Forecast : 0.0189024
Ljung-Box
Test p-value for Missing Person - ARIMA Forecast : 0.3000926
Ljung-Box
Test p-value for STL Residuals: Missing Person : 0.000183457 Processing
category: Motor Vehicle Theft
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1703.658 ARIMA(0,1,0) with drift : 1805.419 ARIMA(1,1,0)(1,0,0)[12] with drift : 1736.452 ARIMA(0,1,1)(0,0,1)[12] with drift : 1708.727 ARIMA(0,1,0) : 1803.394 ARIMA(2,1,2)(0,0,1)[12] with drift : 1711.618 ARIMA(2,1,2)(1,0,0)[12] with drift : 1705.231 ARIMA(2,1,2)(2,0,1)[12] with drift : 1703.142 ARIMA(2,1,2)(2,0,0)[12] with drift : 1702.686 ARIMA(1,1,2)(2,0,0)[12] with drift : 1703.437 ARIMA(2,1,1)(2,0,0)[12] with drift : 1700.583 ARIMA(2,1,1)(1,0,0)[12] with drift : 1704.421 ARIMA(2,1,1)(2,0,1)[12] with drift : 1701.033 ARIMA(2,1,1)(1,0,1)[12] with drift : 1704.269 ARIMA(1,1,1)(2,0,0)[12] with drift : 1701.948 ARIMA(2,1,0)(2,0,0)[12] with drift : 1709.634 ARIMA(3,1,1)(2,0,0)[12] with drift : 1698.399 ARIMA(3,1,1)(1,0,0)[12] with drift : 1709.438 ARIMA(3,1,1)(2,0,1)[12] with drift : 1700.494 ARIMA(3,1,1)(1,0,1)[12] with drift : 1706.779 ARIMA(3,1,0)(2,0,0)[12] with drift : 1709.962 ARIMA(4,1,1)(2,0,0)[12] with drift : 1703.331 ARIMA(3,1,2)(2,0,0)[12] with drift : 1699.469 ARIMA(4,1,0)(2,0,0)[12] with drift : 1711.704 ARIMA(4,1,2)(2,0,0)[12] with drift : 1703.954 ARIMA(3,1,1)(2,0,0)[12] : 1696.935 ARIMA(3,1,1)(1,0,0)[12] : 1707.343 ARIMA(3,1,1)(2,0,1)[12] : 1699.035 ARIMA(3,1,1)(1,0,1)[12] : 1704.673 ARIMA(2,1,1)(2,0,0)[12] : 1698.604 ARIMA(3,1,0)(2,0,0)[12] : 1707.893 ARIMA(4,1,1)(2,0,0)[12] : 1701.831 ARIMA(3,1,2)(2,0,0)[12] : 1697.918 ARIMA(2,1,0)(2,0,0)[12] : 1707.565 ARIMA(2,1,2)(2,0,0)[12] : 1700.685 ARIMA(4,1,0)(2,0,0)[12] : 1709.644 ARIMA(4,1,2)(2,0,0)[12] : 1702.471
Now re-fitting the best model(s) without approximations…
ARIMA(3,1,1)(2,0,0)[12] : 1709.192
Best model: ARIMA(3,1,1)(2,0,0)[12]
Ljung-Box
Test p-value for Motor Vehicle Theft - Naive Forecast : 5.401332e-07
Ljung-Box
Test p-value for Motor Vehicle Theft - Mean Forecast : 0
Ljung-Box
Test p-value for Motor Vehicle Theft - Holt-Winters Forecast : 0.2059866
Ljung-Box
Test p-value for Motor Vehicle Theft - ARIMA Forecast : 0.9971646
Ljung-Box
Test p-value for STL Residuals: Motor Vehicle Theft : 0.06283549
Processing category:
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1775.655 ARIMA(0,1,0) with drift : 1847.328 ARIMA(1,1,0)(1,0,0)[12] with drift : 1840.079 ARIMA(0,1,1)(0,0,1)[12] with drift : 1797.562 ARIMA(0,1,0) : 1845.302 ARIMA(2,1,2)(0,0,1)[12] with drift : 1766.975 ARIMA(2,1,2) with drift : 1764.925 ARIMA(2,1,2)(1,0,0)[12] with drift : 1773.557 ARIMA(1,1,2) with drift : 1766.05 ARIMA(2,1,1) with drift : 1767.775 ARIMA(3,1,2) with drift : Inf ARIMA(2,1,3) with drift : 1766.966 ARIMA(1,1,1) with drift : 1767.053 ARIMA(1,1,3) with drift : 1768.135 ARIMA(3,1,1) with drift : 1768.799 ARIMA(3,1,3) with drift : Inf ARIMA(2,1,2) : 1764.077 ARIMA(2,1,2)(1,0,0)[12] : 1772.894 ARIMA(2,1,2)(0,0,1)[12] : 1766.121 ARIMA(2,1,2)(1,0,1)[12] : 1774.956 ARIMA(1,1,2) : 1764.978 ARIMA(2,1,1) : 1766.719 ARIMA(3,1,2) : 1766.845 ARIMA(2,1,3) : 1766.025 ARIMA(1,1,1) : 1766.393 ARIMA(1,1,3) : 1767.006 ARIMA(3,1,1) : 1768.048 ARIMA(3,1,3) : 1770.041
Now re-fitting the best model(s) without approximations…
ARIMA(2,1,2) : 1767.102
Best model: ARIMA(2,1,2)
Ljung-Box
Test p-value for - Naive Forecast : 9.617126e-05
Ljung-Box
Test p-value for - Mean Forecast : 0
Ljung-Box
Test p-value for - Holt-Winters Forecast : 0.1314561
Ljung-Box
Test p-value for - ARIMA Forecast : 0.8698869
Ljung-Box
Test p-value for STL Residuals: : 1.265572e-05 Processing category:
Weapons
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1003.524 ARIMA(0,1,0) with drift : 1162.605 ARIMA(1,1,0)(1,0,0)[12] with drift : 1092.086 ARIMA(0,1,1)(0,0,1)[12] with drift : 1025.88 ARIMA(0,1,0) : 1160.588 ARIMA(2,1,2)(0,0,1)[12] with drift : 998.7729 ARIMA(2,1,2) with drift : 996.9786 ARIMA(2,1,2)(1,0,0)[12] with drift : 1001.883 ARIMA(1,1,2) with drift : 994.5858 ARIMA(1,1,2)(1,0,0)[12] with drift : 1002.848 ARIMA(1,1,2)(0,0,1)[12] with drift : 996.2263 ARIMA(1,1,2)(1,0,1)[12] with drift : 1002.23 ARIMA(0,1,2) with drift : 1025.601 ARIMA(1,1,1) with drift : 992.578 ARIMA(1,1,1)(1,0,0)[12] with drift : 1000.902 ARIMA(1,1,1)(0,0,1)[12] with drift : 994.2099 ARIMA(1,1,1)(1,0,1)[12] with drift : 1002.827 ARIMA(0,1,1) with drift : 1024.115 ARIMA(1,1,0) with drift : 1087.093 ARIMA(2,1,1) with drift : 1001.741 ARIMA(2,1,0) with drift : 1048.814 ARIMA(1,1,1) : 992.7667
Now re-fitting the best model(s) without approximations…
ARIMA(1,1,1) with drift : 1004.528
Best model: ARIMA(1,1,1) with drift
Ljung-Box
Test p-value for Weapons - Naive Forecast : 1.386472e-08
Ljung-Box
Test p-value for Weapons - Mean Forecast : 0
Ljung-Box
Test p-value for Weapons - Holt-Winters Forecast : 0.9054433
Ljung-Box
Test p-value for Weapons - ARIMA Forecast : 0.9017289 Ljung-Box Test
p-value for STL Residuals: Weapons : 0.05359485 Processing category:
Traffic Except DWI
## Warning in HoltWinters(ts_category): optimization difficulties: ERROR:
## ABNORMAL_TERMINATION_IN_LNSRCH
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1580.283 ARIMA(0,1,0) with drift : 1667.695 ARIMA(1,1,0)(1,0,0)[12] with drift : 1612.382 ARIMA(0,1,1)(0,0,1)[12] with drift : 1563.861 ARIMA(0,1,0) : 1665.675 ARIMA(0,1,1) with drift : 1564.065 ARIMA(0,1,1)(1,0,1)[12] with drift : 1576.074 ARIMA(0,1,1)(0,0,2)[12] with drift : 1564.506 ARIMA(0,1,1)(1,0,0)[12] with drift : 1574.247 ARIMA(0,1,1)(1,0,2)[12] with drift : 1576.949 ARIMA(0,1,0)(0,0,1)[12] with drift : 1668.557 ARIMA(1,1,1)(0,0,1)[12] with drift : 1566.512 ARIMA(0,1,2)(0,0,1)[12] with drift : 1565.105 ARIMA(1,1,0)(0,0,1)[12] with drift : 1602.224 ARIMA(1,1,2)(0,0,1)[12] with drift : 1568.282 ARIMA(0,1,1)(0,0,1)[12] : 1562.327 ARIMA(0,1,1) : 1562.457 ARIMA(0,1,1)(1,0,1)[12] : 1574.509 ARIMA(0,1,1)(0,0,2)[12] : 1562.93 ARIMA(0,1,1)(1,0,0)[12] : 1572.679 ARIMA(0,1,1)(1,0,2)[12] : 1575.276 ARIMA(0,1,0)(0,0,1)[12] : 1666.524 ARIMA(1,1,1)(0,0,1)[12] : 1564.987 ARIMA(0,1,2)(0,0,1)[12] : 1563.52 ARIMA(1,1,0)(0,0,1)[12] : 1600.251 ARIMA(1,1,2)(0,0,1)[12] : 1566.702
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1)(0,0,1)[12] : 1566.002
Best model: ARIMA(0,1,1)(0,0,1)[12]
Ljung-Box
Test p-value for Traffic Except DWI - Naive Forecast : 6.126916e-09
Ljung-Box
Test p-value for Traffic Except DWI - Mean Forecast : 0
Ljung-Box
Test p-value for Traffic Except DWI - Holt-Winters Forecast : 0.73422
Ljung-Box
Test p-value for Traffic Except DWI - ARIMA Forecast : 0.8634272
Ljung-Box
Test p-value for STL Residuals: Traffic Except DWI : 0.1191584
Processing category: Aggravated Assault
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1177.156 ARIMA(0,1,0) with drift : 1361.828 ARIMA(1,1,0)(1,0,0)[12] with drift : 1286.102 ARIMA(0,1,1)(0,0,1)[12] with drift : 1166.586 ARIMA(0,1,0) : 1359.802 ARIMA(0,1,1) with drift : 1165.956 ARIMA(0,1,1)(1,0,0)[12] with drift : 1185.043 ARIMA(0,1,1)(1,0,1)[12] with drift : 1183.359 ARIMA(1,1,1) with drift : 1169.105 ARIMA(0,1,2) with drift : 1167.735 ARIMA(1,1,0) with drift : 1279.598 ARIMA(1,1,2) with drift : 1169.528 ARIMA(0,1,1) : 1164.803 ARIMA(0,1,1)(1,0,0)[12] : 1183.678 ARIMA(0,1,1)(0,0,1)[12] : 1165.373 ARIMA(0,1,1)(1,0,1)[12] : 1182.282 ARIMA(1,1,1) : 1168.007 ARIMA(0,1,2) : 1166.431 ARIMA(1,1,0) : 1277.558 ARIMA(1,1,2) : 1168.936
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1) : 1167.382
Best model: ARIMA(0,1,1)
Ljung-Box
Test p-value for Aggravated Assault - Naive Forecast : 1.453305e-10
Ljung-Box
Test p-value for Aggravated Assault - Mean Forecast : 0.04504411
Ljung-Box
Test p-value for Aggravated Assault - Holt-Winters Forecast : 0.773886
Ljung-Box
Test p-value for Aggravated Assault - ARIMA Forecast : 0.7395607
Ljung-Box
Test p-value for STL Residuals: Aggravated Assault : 0.01247429
Processing category: Stolen Property
Fitting models using approximations to speed things up…
ARIMA(2,0,2)(1,0,1)[12] with non-zero mean : 714.1806 ARIMA(0,0,0) with non-zero mean : 711.3391 ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 711.3097 ARIMA(0,0,1)(0,0,1)[12] with non-zero mean : 714.7307 ARIMA(0,0,0) with zero mean : 1027.13 ARIMA(1,0,0) with non-zero mean : 713.6401 ARIMA(1,0,0)(2,0,0)[12] with non-zero mean : 711.7651 ARIMA(1,0,0)(1,0,1)[12] with non-zero mean : 712.8517 ARIMA(1,0,0)(0,0,1)[12] with non-zero mean : 715.0504 ARIMA(1,0,0)(2,0,1)[12] with non-zero mean : 712.3669 ARIMA(0,0,0)(1,0,0)[12] with non-zero mean : 708.9065 ARIMA(0,0,0)(2,0,0)[12] with non-zero mean : 709.1305 ARIMA(0,0,0)(1,0,1)[12] with non-zero mean : 710.4019 ARIMA(0,0,0)(0,0,1)[12] with non-zero mean : 712.6614 ARIMA(0,0,0)(2,0,1)[12] with non-zero mean : 710.1593 ARIMA(0,0,1)(1,0,0)[12] with non-zero mean : 710.9787 ARIMA(1,0,1)(1,0,0)[12] with non-zero mean : 712.4647 ARIMA(0,0,0)(1,0,0)[12] with zero mean : 826.1549
Now re-fitting the best model(s) without approximations…
ARIMA(0,0,0)(1,0,0)[12] with non-zero mean : 712.5003
Best model: ARIMA(0,0,0)(1,0,0)[12] with non-zero mean
Ljung-Box
Test p-value for Stolen Property - Naive Forecast : 5.844937e-08
Ljung-Box
Test p-value for Stolen Property - Mean Forecast : 0.4597357
Ljung-Box
Test p-value for Stolen Property - Holt-Winters Forecast : 0.7589937
Ljung-Box
Test p-value for Stolen Property - ARIMA Forecast : 0.5038815
Ljung-Box
Test p-value for STL Residuals: Stolen Property : 0.09713705 Processing
category: Alcohol Offenses
Fitting models using approximations to speed things up…
ARIMA(2,0,2)(1,0,1)[12] with non-zero mean : 627.0149 ARIMA(0,0,0) with non-zero mean : 650.9832 ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 626.2375 ARIMA(0,0,1)(0,0,1)[12] with non-zero mean : 649.8751 ARIMA(0,0,0) with zero mean : 901.3115 ARIMA(1,0,0) with non-zero mean : 647.0568 ARIMA(1,0,0)(2,0,0)[12] with non-zero mean : 635.2371 ARIMA(1,0,0)(1,0,1)[12] with non-zero mean : 628.0333 ARIMA(1,0,0)(0,0,1)[12] with non-zero mean : 648.9856 ARIMA(1,0,0)(2,0,1)[12] with non-zero mean : 635.9113 ARIMA(0,0,0)(1,0,0)[12] with non-zero mean : 632.6926 ARIMA(2,0,0)(1,0,0)[12] with non-zero mean : 623.3495 ARIMA(2,0,0) with non-zero mean : 642.599 ARIMA(2,0,0)(2,0,0)[12] with non-zero mean : 632.8493 ARIMA(2,0,0)(1,0,1)[12] with non-zero mean : 624.7798 ARIMA(2,0,0)(0,0,1)[12] with non-zero mean : 644.6912 ARIMA(2,0,0)(2,0,1)[12] with non-zero mean : 634.7545 ARIMA(3,0,0)(1,0,0)[12] with non-zero mean : 624.2646 ARIMA(2,0,1)(1,0,0)[12] with non-zero mean : 623.4951 ARIMA(1,0,1)(1,0,0)[12] with non-zero mean : 620.6952 ARIMA(1,0,1) with non-zero mean : 641.1887 ARIMA(1,0,1)(2,0,0)[12] with non-zero mean : 629.1078 ARIMA(1,0,1)(1,0,1)[12] with non-zero mean : 621.9377 ARIMA(1,0,1)(0,0,1)[12] with non-zero mean : 643.2829 ARIMA(1,0,1)(2,0,1)[12] with non-zero mean : 631.2013 ARIMA(0,0,1)(1,0,0)[12] with non-zero mean : 628.2296 ARIMA(1,0,2)(1,0,0)[12] with non-zero mean : 622.1069 ARIMA(0,0,2)(1,0,0)[12] with non-zero mean : 625.8714 ARIMA(2,0,2)(1,0,0)[12] with non-zero mean : 625.6268 ARIMA(1,0,1)(1,0,0)[12] with zero mean : Inf
Now re-fitting the best model(s) without approximations…
ARIMA(1,0,1)(1,0,0)[12] with non-zero mean : 642.3412
Best model: ARIMA(1,0,1)(1,0,0)[12] with non-zero mean
Ljung-Box
Test p-value for Alcohol Offenses - Naive Forecast : 7.392966e-07
Ljung-Box
Test p-value for Alcohol Offenses - Mean Forecast : 0.02706037
Ljung-Box
Test p-value for Alcohol Offenses - Holt-Winters Forecast : 0.6155893
Ljung-Box
Test p-value for Alcohol Offenses - ARIMA Forecast : 0.9316562
Ljung-Box
Test p-value for STL Residuals: Alcohol Offenses : 0.2785603 Processing
category: Drunkenness
Fitting models using approximations to speed things up…
ARIMA(2,0,2)(1,0,1)[12] with non-zero mean : 685.354 ARIMA(0,0,0) with non-zero mean : 679.524 ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 680.934 ARIMA(0,0,1)(0,0,1)[12] with non-zero mean : 674.5013 ARIMA(0,0,0) with zero mean : 980.1668 ARIMA(0,0,1) with non-zero mean : 674.0378 ARIMA(0,0,1)(1,0,0)[12] with non-zero mean : 681.7818 ARIMA(0,0,1)(1,0,1)[12] with non-zero mean : 683.6758 ARIMA(1,0,1) with non-zero mean : 670.7251 ARIMA(1,0,1)(1,0,0)[12] with non-zero mean : 678.9194 ARIMA(1,0,1)(0,0,1)[12] with non-zero mean : 670.2638 ARIMA(1,0,1)(1,0,1)[12] with non-zero mean : 680.0316 ARIMA(1,0,1)(0,0,2)[12] with non-zero mean : 672.1654 ARIMA(1,0,1)(1,0,2)[12] with non-zero mean : 681.9477 ARIMA(1,0,0)(0,0,1)[12] with non-zero mean : 673.4422 ARIMA(2,0,1)(0,0,1)[12] with non-zero mean : 673.0538 ARIMA(1,0,2)(0,0,1)[12] with non-zero mean : 672.1487 ARIMA(0,0,0)(0,0,1)[12] with non-zero mean : 679.2694 ARIMA(0,0,2)(0,0,1)[12] with non-zero mean : 673.9827 ARIMA(2,0,0)(0,0,1)[12] with non-zero mean : 674.22 ARIMA(2,0,2)(0,0,1)[12] with non-zero mean : 675.1838 ARIMA(1,0,1)(0,0,1)[12] with zero mean : Inf
Now re-fitting the best model(s) without approximations…
ARIMA(1,0,1)(0,0,1)[12] with non-zero mean : 671.1076
Best model: ARIMA(1,0,1)(0,0,1)[12] with non-zero mean
Ljung-Box
Test p-value for Drunkenness - Naive Forecast : 2.781361e-08
Ljung-Box
Test p-value for Drunkenness - Mean Forecast : 0.00274216
Ljung-Box
Test p-value for Drunkenness - Holt-Winters Forecast : 0.04449286
Ljung-Box
Test p-value for Drunkenness - ARIMA Forecast : 0.4497937
Ljung-Box
Test p-value for STL Residuals: Drunkenness : 0.0002278958 Processing
category: DUI
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1695.421 ARIMA(0,1,0) with drift : 1761.531 ARIMA(1,1,0)(1,0,0)[12] with drift : 1706.67 ARIMA(0,1,1)(0,0,1)[12] with drift : 1675.868 ARIMA(0,1,0) : 1759.504 ARIMA(0,1,1) with drift : 1678.236 ARIMA(0,1,1)(1,0,1)[12] with drift : 1688.819 ARIMA(0,1,1)(0,0,2)[12] with drift : 1676.91 ARIMA(0,1,1)(1,0,0)[12] with drift : 1686.84 ARIMA(0,1,1)(1,0,2)[12] with drift : 1688.73 ARIMA(0,1,0)(0,0,1)[12] with drift : 1757.175 ARIMA(1,1,1)(0,0,1)[12] with drift : 1677.86 ARIMA(0,1,2)(0,0,1)[12] with drift : 1677.427 ARIMA(1,1,0)(0,0,1)[12] with drift : 1696.446 ARIMA(1,1,2)(0,0,1)[12] with drift : 1679.656 ARIMA(0,1,1)(0,0,1)[12] : 1673.847 ARIMA(0,1,1) : 1676.258 ARIMA(0,1,1)(1,0,1)[12] : 1686.794 ARIMA(0,1,1)(0,0,2)[12] : 1674.866 ARIMA(0,1,1)(1,0,0)[12] : 1684.831 ARIMA(0,1,1)(1,0,2)[12] : 1686.712 ARIMA(0,1,0)(0,0,1)[12] : 1755.132 ARIMA(1,1,1)(0,0,1)[12] : 1675.841 ARIMA(0,1,2)(0,0,1)[12] : 1675.388 ARIMA(1,1,0)(0,0,1)[12] : 1694.398 ARIMA(1,1,2)(0,0,1)[12] : 1677.63
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1)(0,0,1)[12] : 1677.603
Best model: ARIMA(0,1,1)(0,0,1)[12]
Ljung-Box
Test p-value for DUI - Naive Forecast : 3.528362e-10
Ljung-Box
Test p-value for DUI - Mean Forecast : 0
Ljung-Box
Test p-value for DUI - Holt-Winters Forecast : 0.1799627
Ljung-Box
Test p-value for DUI - ARIMA Forecast : 0.2390681
Ljung-Box
Test p-value for STL Residuals: DUI : 0.002544789 Processing category:
Suicide
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : 1291.187 ARIMA(0,1,0) with drift : 1459.012 ARIMA(1,1,0)(1,0,0)[12] with drift : 1384.184 ARIMA(0,1,1)(0,0,1)[12] with drift : 1278.22 ARIMA(0,1,0) : 1456.988 ARIMA(0,1,1) with drift : 1276.825 ARIMA(0,1,1)(1,0,0)[12] with drift : 1279.729 ARIMA(0,1,1)(1,0,1)[12] with drift : 1279.968 ARIMA(1,1,1) with drift : 1282.467 ARIMA(0,1,2) with drift : 1278.88 ARIMA(1,1,0) with drift : 1378.819 ARIMA(1,1,2) with drift : 1278.952 ARIMA(0,1,1) : 1275.03 ARIMA(0,1,1)(1,0,0)[12] : 1277.722 ARIMA(0,1,1)(0,0,1)[12] : 1276.411 ARIMA(0,1,1)(1,0,1)[12] : 1277.909 ARIMA(1,1,1) : 1280.433 ARIMA(0,1,2) : 1277.072 ARIMA(1,1,0) : 1376.779 ARIMA(1,1,2) : 1276.893
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1) : 1276.485
Best model: ARIMA(0,1,1)
Ljung-Box
Test p-value for Suicide - Naive Forecast : 3.306901e-10
Ljung-Box
Test p-value for Suicide - Mean Forecast : 0.0008676618
Ljung-Box
Test p-value for Suicide - Holt-Winters Forecast : 0.9879291
Ljung-Box
Test p-value for Suicide - ARIMA Forecast : 0.9654798
Ljung-Box
Test p-value for STL Residuals: Suicide : 0.1157256 Processing category:
Arson
ARIMA(2,0,2)(1,0,1)[12] with non-zero mean : 281.9897 ARIMA(0,0,0) with non-zero mean : 272.8251 ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 276.4537 ARIMA(0,0,1)(0,0,1)[12] with non-zero mean : 276.3587 ARIMA(0,0,0) with zero mean : 446.7236 ARIMA(0,0,0)(1,0,0)[12] with non-zero mean : 274.305 ARIMA(0,0,0)(0,0,1)[12] with non-zero mean : 274.2119 ARIMA(0,0,0)(1,0,1)[12] with non-zero mean : 276.2835 ARIMA(1,0,0) with non-zero mean : 274.9216 ARIMA(0,0,1) with non-zero mean : 274.9218 ARIMA(1,0,1) with non-zero mean : 275.7804
Best model: ARIMA(0,0,0) with non-zero mean
Ljung-Box
Test p-value for Arson - Naive Forecast : 0.0001476746
Ljung-Box
Test p-value for Arson - Mean Forecast : 0.2801111
Ljung-Box
Test p-value for Arson - Holt-Winters Forecast : 0.5006736
Ljung-Box
Test p-value for Arson - ARIMA Forecast : 0.2801111
Ljung-Box
Test p-value for STL Residuals: Arson : 0.05347482 Processing category:
Embezzlement
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : Inf ARIMA(0,1,0) with drift : 1256.522 ARIMA(1,1,0)(1,0,0)[12] with drift : 1166.802 ARIMA(0,1,1)(0,0,1)[12] with drift : 1083.445 ARIMA(0,1,0) : 1254.494 ARIMA(0,1,1) with drift : 1081.682 ARIMA(0,1,1)(1,0,0)[12] with drift : 1065.678 ARIMA(0,1,1)(2,0,0)[12] with drift : 1060.513 ARIMA(0,1,1)(2,0,1)[12] with drift : Inf ARIMA(0,1,1)(1,0,1)[12] with drift : Inf ARIMA(0,1,0)(2,0,0)[12] with drift : 1247.022 ARIMA(1,1,1)(2,0,0)[12] with drift : 1064.631 ARIMA(0,1,2)(2,0,0)[12] with drift : 1060.997 ARIMA(1,1,0)(2,0,0)[12] with drift : 1166.883 ARIMA(1,1,2)(2,0,0)[12] with drift : 1066.862 ARIMA(0,1,1)(2,0,0)[12] : 1061.268
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1)(2,0,0)[12] with drift : Inf ARIMA(0,1,2)(2,0,0)[12] with drift : Inf ARIMA(0,1,1)(2,0,0)[12] : 1077.052
Best model: ARIMA(0,1,1)(2,0,0)[12]
Ljung-Box
Test p-value for Embezzlement - Naive Forecast : 1.563849e-11
Ljung-Box
Test p-value for Embezzlement - Mean Forecast : 0
Ljung-Box
Test p-value for Embezzlement - Holt-Winters Forecast : 0.7276434
Ljung-Box
Test p-value for Embezzlement - ARIMA Forecast : 0.5963342
Ljung-Box
Test p-value for STL Residuals: Embezzlement : 0.02397306 Processing
category: Disorderly Conduct
Fitting models using approximations to speed things up…
ARIMA(2,1,2)(1,0,1)[12] with drift : Inf ARIMA(0,1,0) with drift : 819.5349 ARIMA(1,1,0)(1,0,0)[12] with drift : 766.5201 ARIMA(0,1,1)(0,0,1)[12] with drift : 700.2158 ARIMA(0,1,0) : 817.4989 ARIMA(0,1,1) with drift : 698.1954 ARIMA(0,1,1)(1,0,0)[12] with drift : 705.081 ARIMA(0,1,1)(1,0,1)[12] with drift : 705.1481 ARIMA(1,1,1) with drift : 700.3085 ARIMA(0,1,2) with drift : 699.4031 ARIMA(1,1,0) with drift : 760.4518 ARIMA(1,1,2) with drift : 701.9666 ARIMA(0,1,1) : 696.1994 ARIMA(0,1,1)(1,0,0)[12] : 703.0352 ARIMA(0,1,1)(0,0,1)[12] : 698.2004 ARIMA(0,1,1)(1,0,1)[12] : 703.0692 ARIMA(1,1,1) : 698.2807 ARIMA(0,1,2) : 697.3876 ARIMA(1,1,0) : 758.3982 ARIMA(1,1,2) : 699.9114
Now re-fitting the best model(s) without approximations…
ARIMA(0,1,1) : 695.5634
Best model: ARIMA(0,1,1)
Ljung-Box
Test p-value for Disorderly Conduct - Naive Forecast : 2.631229e-14
Ljung-Box
Test p-value for Disorderly Conduct - Mean Forecast : 0.0005170751
Ljung-Box
Test p-value for Disorderly Conduct - Holt-Winters Forecast : 0.07636961
Ljung-Box
Test p-value for Disorderly Conduct - ARIMA Forecast : 0.06530758
Ljung-Box
Test p-value for STL Residuals: Disorderly Conduct : 0.001545906
Processing category: Pornography
ARIMA(2,0,2)(1,0,1)[12] with non-zero mean : Inf ARIMA(0,0,0) with non-zero mean : 87.43653 ARIMA(1,0,0)(1,0,0)[12] with non-zero mean : 90.56536 ARIMA(0,0,1)(0,0,1)[12] with non-zero mean : 90.5408 ARIMA(0,0,0) with zero mean : 184.6775 ARIMA(0,0,0)(1,0,0)[12] with non-zero mean : 89.59076 ARIMA(0,0,0)(0,0,1)[12] with non-zero mean : 89.62992 ARIMA(0,0,0)(1,0,1)[12] with non-zero mean : Inf ARIMA(1,0,0) with non-zero mean : 88.39065 ARIMA(0,0,1) with non-zero mean : 88.28496 ARIMA(1,0,1) with non-zero mean : Inf
Best model: ARIMA(0,0,0) with non-zero mean
Ljung-Box
Test p-value for Pornography - Naive Forecast : 0.005772429
Ljung-Box
Test p-value for Pornography - Mean Forecast : 0.6229255
Ljung-Box
Test p-value for Pornography - Holt-Winters Forecast : 0.8459962
Ljung-Box
Test p-value for Pornography - ARIMA Forecast : 0.6229255
Ljung-Box
Test p-value for STL Residuals: Pornography : 0.07096858 Processing
category: Murder Insufficient data for STL decomposition for category:
Murder Processing category: Non Negligent Traff Insufficient data for
STL decomposition for category: Non Negligent Traff Processing category:
Gambling Insufficient data for STL decomposition for category: Gambling
Processing category: Calls for Service
ARIMA(2,0,2) with non-zero mean : Inf ARIMA(0,0,0) with non-zero mean : 8.772111 ARIMA(1,0,0) with non-zero mean : 11.13465 ARIMA(0,0,1) with non-zero mean : 11.09062 ARIMA(0,0,0) with zero mean : 84.20079 ARIMA(1,0,1) with non-zero mean : Inf
Best model: ARIMA(0,0,0) with non-zero mean
Ljung-Box
Test p-value for Calls for Service - Naive Forecast : 0.0006768923
Ljung-Box
Test p-value for Calls for Service - Mean Forecast : 0.1919992
Ljung-Box
Test p-value for Calls for Service - Holt-Winters Forecast : NA
Ljung-Box
Test p-value for Calls for Service - ARIMA Forecast : 0.1919992
Ljung-Box
Test p-value for STL Residuals: Calls for Service : 0.004584851
Processing category: Non Criminal Detain
ARIMA(2,1,2) with drift : 199.0386 ARIMA(0,1,0) with drift : 208.4051 ARIMA(1,1,0) with drift : 197.0616 ARIMA(0,1,1) with drift : 196.612 ARIMA(0,1,0) : 206.1629 ARIMA(1,1,1) with drift : 198.1586 ARIMA(0,1,2) with drift : 198.2635 ARIMA(1,1,2) with drift : 196.5745 ARIMA(1,1,3) with drift : 197.847 ARIMA(0,1,3) with drift : 201.0605 ARIMA(2,1,1) with drift : 200.6467 ARIMA(2,1,3) with drift : 198.7708 ARIMA(1,1,2) : 194.3352 ARIMA(0,1,2) : 196.0607 ARIMA(1,1,1) : 196.0295 ARIMA(2,1,2) : 196.7686 ARIMA(1,1,3) : 195.7849 ARIMA(0,1,1) : 194.9304 ARIMA(0,1,3) : 198.3212 ARIMA(2,1,1) : 198.4078 ARIMA(2,1,3) : Inf
Best model: ARIMA(1,1,2)
Ljung-Box
Test p-value for Non Criminal Detain - Naive Forecast : 3.462545e-07
Ljung-Box
Test p-value for Non Criminal Detain - Mean Forecast : 1.180947e-08
Ljung-Box
Test p-value for Non Criminal Detain - Holt-Winters Forecast : 0.3198892
Ljung-Box
Test p-value for Non Criminal Detain - ARIMA Forecast : 0.7445539
Ljung-Box
Test p-value for STL Residuals: Non Criminal Detain : 0.02024312
Processing category: NA Error processing category: NA - no rows to
aggregate Successfully processed categories: Larceny Vandalism Forgery
Simple Assault Burglary Fraud All Other Offenses Drugs Missing Person
Motor Vehicle Theft Weapons Traffic Except DWI Aggravated Assault Stolen
Property Alcohol Offenses Drunkenness DUI Suicide Arson Embezzlement
Disorderly Conduct Pornography Calls for Service Non Criminal Detain
Categories with issues:
RMSE for category: Larceny Naive RMSE: 25.23893 Mean RMSE: 23.41082 Holt-Winters RMSE: 20.45496 STL RMSE: 7.336082 ARIMA RMSE: 21.13999
RMSE for category: Vandalism Naive RMSE: 11.98034 Mean RMSE: 13.5175 Holt-Winters RMSE: 9.69951 STL RMSE: 3.334571 ARIMA RMSE: 9.60029
RMSE for category: Forgery Naive RMSE: 4.284809 Mean RMSE: 5.195624 Holt-Winters RMSE: 4.417228 STL RMSE: 0.5847482 ARIMA RMSE: 3.914415
RMSE for category: Simple Assault Naive RMSE: 8.395229 Mean RMSE: 7.380232 Holt-Winters RMSE: 6.6264 STL RMSE: 1.985979 ARIMA RMSE: 6.365325
RMSE for category: Burglary Naive RMSE: 10.23843 Mean RMSE: 12.20652 Holt-Winters RMSE: 9.02529 STL RMSE: 1.610248 ARIMA RMSE: 7.892439
RMSE for category: Fraud Naive RMSE: 10.05488 Mean RMSE: 10.1406 Holt-Winters RMSE: 8.706708 STL RMSE: 1.635293 ARIMA RMSE: 8.192234
RMSE for category: All Other Offenses Naive RMSE: 12.83543 Mean RMSE: 21.06401 Holt-Winters RMSE: 10.45179 STL RMSE: 3.681927 ARIMA RMSE: 10.49191
RMSE for category: Drugs Naive RMSE: 8.704508 Mean RMSE: 10.54179 Holt-Winters RMSE: 7.788851 STL RMSE: 1.216346 ARIMA RMSE: 7.125262
RMSE for category: Missing Person Naive RMSE: 2.625713 Mean RMSE: 2.252828 Holt-Winters RMSE: 2.18036 STL RMSE: 0.5578251 ARIMA RMSE: 2.057782
RMSE for category: Motor Vehicle Theft Naive RMSE: 5.101007 Mean RMSE: 5.950275 Holt-Winters RMSE: 4.714105 STL RMSE: 0.9927359 ARIMA RMSE: 4.222176
RMSE for category: Weapons Naive RMSE: 2.159665 Mean RMSE: 1.764051 Holt-Winters RMSE: 1.70406 STL RMSE: 0.3049746 ARIMA RMSE: 1.577494
RMSE for category: Traffic Except DWI Naive RMSE: 4.59874 Mean RMSE: 8.062352 Holt-Winters RMSE: 4.046414 STL RMSE: 0.6869745 ARIMA RMSE: 3.797009
RMSE for category: Aggravated Assault Naive RMSE: 2.502233 Mean RMSE: 1.799438 Holt-Winters RMSE: 1.91639 STL RMSE: 0.1886339 ARIMA RMSE: 1.775804
RMSE for category: Stolen Property Naive RMSE: 1.617634 Mean RMSE: 1.141253 Holt-Winters RMSE: 1.284778 STL RMSE: 0.2934632 ARIMA RMSE: 1.138896
RMSE for category: Alcohol Offenses Naive RMSE: 1.571527 Mean RMSE: 1.22936 Holt-Winters RMSE: 1.392433 STL RMSE: 0.2721775 ARIMA RMSE: 1.183851
RMSE for category: Drunkenness Naive RMSE: 1.36626 Mean RMSE: 1.078385 Holt-Winters RMSE: 1.16832 STL RMSE: 0.2647226 ARIMA RMSE: 1.043563
RMSE for category: DUI Naive RMSE: 5.371744 Mean RMSE: 10.79571 Holt-Winters RMSE: 4.745643 STL RMSE: 0.9076895 ARIMA RMSE: 4.577422
RMSE for category: Suicide Naive RMSE: 2.883804 Mean RMSE: 2.147559 Holt-Winters RMSE: 2.208467 STL RMSE: 0.3338795 ARIMA RMSE: 2.095057
RMSE for category: Arson Naive RMSE: 1.101086 Mean RMSE: 0.7863486 Holt-Winters RMSE: 0.8619329 STL RMSE: 0.2284818 ARIMA RMSE: 0.7863486
RMSE for category: Embezzlement Naive RMSE: 2.580039 Mean RMSE: 2.053076 Holt-Winters RMSE: 1.989035 STL RMSE: 0.4795919 ARIMA RMSE: 1.806899
RMSE for category: Disorderly Conduct Naive RMSE: 1.472461 Mean RMSE: 1.130143 Holt-Winters RMSE: 1.235971 STL RMSE: 0.1677278 ARIMA RMSE: 1.109612
RMSE for category: Pornography Naive RMSE: 0.8009428 Mean RMSE: 0.5227998 Holt-Winters RMSE: 0.6276156 STL RMSE: 0.3017362 ARIMA RMSE: 0.5227998
RMSE for category: Calls for Service Naive RMSE: 0.3922323 Mean RMSE: 0.2618914 Holt-Winters RMSE: 0.2588865 STL RMSE: 0.1847502 ARIMA RMSE: 0.2618914
RMSE for category: Non Criminal Detain Naive RMSE: 4.862824 Mean RMSE: 5.339571 Holt-Winters RMSE: 4.50186 STL RMSE: 2.273436 ARIMA RMSE: 3.577469 Overall, the forecast suggests a slight decline in crime rates in the near future, followed by minor fluctuations. The trend indicates that crime levels are likely to stabilize and remain close to the current levels by the end of the year. r The RMSE values from different forecasting methods highlight their predictive accuracy for the given time series data. Among the methods, the Moving Average with a window size of 3 (MA3) achieves the lowest RMSE (28.70), suggesting it provides the most accurate forecasts. The Holt-Winters (HW) method also performs well (RMSE: 42.97), capturing trends and seasonality effectively. In contrast, the Naive method has the highest RMSE (55.76), indicating less reliable predictions. These results emphasize the importance of selecting the most suitable method, like MA3, for forecasting in this context. Based on the individual forecasts for various crime categories, the analysis reveals distinct trends. Larceny, forgery, drugs, missing person, motor vehicle theft, and stolen property initially show a decrease, with stolen property subsequently increasing. Vandalism, burglary, fraud, weapons offenses, and suicide are projected to rise, albeit at varying rates. Simple assault, alcohol-related offenses, DUI, aggravated assault, arson, embezzlement, and other offenses are expected to remain stable, with minor fluctuations in some cases. Traffic violations, excluding DWI, are predicted to either decline slightly or remain steady. These insights provide a nuanced understanding of potential future trends in crime rates across different categories. The RMSE values highlight the accuracy of different forecasting methods, with lower values indicating better performance. STL decomposition consistently achieves the lowest RMSE across most categories, making it the most reliable method, particularly for crimes like Larceny, Vandalism, and Forgery. ARIMA and Holt-Winters also perform well in capturing trends and seasonality, while Naive and Mean forecasts generally have higher RMSE, indicating lower accuracy. These comparisons guide the selection of effective models for forecasting crime trends.
Call: lm(formula = Lon ~ Lat, data = df)
Residuals: Min 1Q Median 3Q Max -0.274323 -0.029408 0.008839 0.032409 0.135637
Coefficients: Estimate Std. Error t value Pr(>|t|)
(Intercept) -56.089368 0.124733 -449.7 <2e-16 Lat
-0.634565 0.003486 -182.0 <2e-16 — Signif. codes: 0
‘’ 0.001 ’’ 0.01 ’’ 0.05 ‘.’ 0.1 ’ ’ 1
Residual standard error: 0.03845 on 115090 degrees of freedom (6024 observations deleted due to missingness) Multiple R-squared: 0.2235, Adjusted R-squared: 0.2235 F-statistic: 3.313e+04 on 1 and 115090 DF, p-value: < 2.2e-16
The regression model shows a significant relationship between latitude and longitude, with an R-squared of 22.35%, indicating that 22.35% of longitude variance is explained by latitude. The coefficients are highly significant (p < 2.2e-16), and the residual standard error is low (0.03845), suggesting a strong but not comprehensive fit.
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 6024 rows containing non-finite outside the scale range
## (`stat_smooth()`).
## Warning: Removed 6024 rows containing missing values or values outside the scale range
## (`geom_point()`).
Here is performing simple regression on lattitude and longitude, since
the dataset has a lot of categoriacal data, regression might not be a
best strategy, however the linear regression model predicts the
longitude (Lon) based on the latitude (Lat) of crime incidents. The
scatter plot visualizes the relationship between Lat and Lon, with the
regression line (in red) showing the best fit. The summary of the
regression model includes coefficients, their significance, and overall
model metrics, providing insights into the spatial distribution of
crimes.